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A Josephson junction consisting of an Anderson impurity weakly coupled with two-band and 
single-band superconductors exposes a time reversal breaking ground state when the coupling be- 
tween the two bands exceeds a certain threshold. The critical regime occurs around local moment 
formation. This indicates a fundamental and distinct role of strong correlations: Driving a system 
into a time reversal breaking ground state. One of the observable consequences is that the impurity 
magnetization in this phase is reduced. 

Motivation: Recently, there is a renewed interest in the physics of multi-band superconductors which are believed to 
encode the physics behind the new Iron based superconductorsi*^!^. Historically, interest in two-band superconductors 
is related to experiments suggesting that in some transition metals, there should be a second mechanism responsible 
for pairing beyond that of electron-phonon interaction^i^ii. For example, in NbaSn the isotopic effect is found to be 
very weak^. Inter-band coupling also provides a natural mechanism for augmenting Tc^^ 

The physics associated with a two-band superconductor can be examined in a Josephson junction involving two-band 
superconductors on one of its sides and a third superconductor on the other sideiS. Within a Ginzburg-Landau 
formalism the order parameters '4'i,il'2 of the two bands are coupled among themselves by a term —JIle[il)iip2]i ^^nd 
to the single-band superconductor (on the opposite side of the junction) with tunneling strengths Fi and F2. A 
time-reversal violating ground state (TRVGS) was found for Fi w F2, due to frustration between tire three order 
parameters, while otherwise, a time reversal conserving ground state (TRCGS) prevails. A natural question which 
motivates the present research is what would happen if tire Josephson junction contains an Anderson impurity. Indeed, 
strong correlations are expected to affect these findings and alter the pertinent physics in a profound way. 
In this work we study a system consisting of an Anderson impurity (level energy e and correlation strength U) 
weakly coupled on the left with a two-band superconductor (order parameters ipn = A„e'^", n = 1,2 and inter-band 
Josepshon coupling J) and on the right to a single-band superconductor (order parameter '03 = Aae**^). This is briefly 
denoted as 2B ~ U — IB Josephson junction (see Fig[T]). A TVRGS occurs if the free energy J^[0 = {61, 62, 0^)] has 




FIG. 1: Geometry of the 2B - U - IB junction 

a minimum at some point 9 for which there is a non-zero Josephson current loop Jc even in the absence of magnetic 
field (see RefJ^ for the case of 2B — IB Josephson junction without Anderson impurity). On the technical part, the 
mean-field (Hartree-Fock) approach to the IB — U — IB problemi^ii^ can be used to treat the 2B — U — IB system 
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as well. It enables the elucidation of the ground-state configuration 0(J, U) together with the free energy !F{0), the 
Josephson current J^dO) the impurity occupation n\jnp[9) and magnetization m^mpid). As it turn out, however, the 
physics is distinct since there are now three order parameters and two of them are coupled. 

Our main result pertains to the nature of the ground state TR symmetry in the U ~ J plane, for high asymmetry 
parameter, (5 = ^ 1, where no TRVGS exists^^ on the line [/ = 0. Two phases are identified, one for which 
there is TRCGS and one for which there is TRVGS, separated by a sharp border line. The TRVGS occurs for 
Ui{5) <U < U2{5) and for J > Jc{U,5) > 0. Here Ui < Uq and U2 > Uq where Uq defined by e + ^ = is the 
particle hole symmetric point and Jc{U,S) is some critical threshold. 
Formalism: The mean-field Hamiltonian is written as, 

H = H^ + HR + Hi + Hi^p. (1) 



The first two terms Hl, Hji describe (within the BCS formalism) the two band superconductor on the left (Hl) and the 
single band superconductor oi 
and Hfj = H3 are written as. 



single band superconductor on the right {Hji). Explicitly, in terms of quasi-particle field operators, Hj^ — X]n=i 



H,,= j ¥^{Y)nn^n{v)dv, ^'„(r) 



(2) 



The two-band Hamiltonian densities 'Hn=i,2 are derived rom the two-band model Hamiltonian of Refi^. In this 
Hamiltonian (equation (1) therein), there is only intraband pairing (and no intcrband pairing) and the electrons in 
the two bands are coupled only through an intcrband Josephson effect (see below). The single band Hamiltonian 
density 7^3 has the standard Bogoliubov-De Gennes structure. Thus, 

\^A„e-*«" J™e-*«- -£„(-zV) + m/' ' \ Age-^^^ ~e3{-iV) + ^i) ' 

In Eqs. ^ m ^ n ~ 1,2 and the various quantities are defined as follows: e(— iV) is the kinetic energy operator derived 
from the corresonding energy dispersion functions £n(k), and fi is the chemical potential. Moreover, A„e'^" = Vnipn 
where ?/;„ = ^ dv{'il^ni{'^)'4'n-\{^)) is the order parameter of superconductor n ~ 1,2,3 and Vn is the corresponding 
strength of the pairing potential. Similarly, Jme*^"" = /i/'m encodes the pairing field in band n ^ m due to electrons 
pairing in band m (intcrband Josephson effect) and / is the strength of the coupling between the two bands^. 
The tunneling part Hi contains hopping between the impurity to each one of the three superconductors (with different 
strengths tn=i,2.3), which occurs at a single point. Finally, the strong correlation part has the usual structure of an 
Anderson impurity Hamiltonian. Explicitly, 



= - E *n[*U0)T3C + CtrgVl/iO) + h.c], C 



H,^^ = eC^T^C+hj[C^C]\ ^=e+|. (4) 



The procedure for calculating the free energy T[6\, 62, 63) is a modified version of the algorithm used in RefsJ^ii^. 
It involves an Eucledean path integral for the partition function Z = e~^'^ in terms of Grassman fields and employing 
Hubbard- Stratonovich transformation for treating the quartic term in i/imp at the expense of an additional integration 
on a new field 7. The latter is carried out within the saddle point approximation leading to a self consistent equation 
(Hartree-Fock approximation). Notice that unlike RefJ^ where the pairing field away from the impurity is also solved 
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self-consistcntly within a Ginzburg-Landau approximation, the pairing fields A„'s are assumed to take a constant 
value in the superconductors. 

Beside the density of states at the Fermi energy -/V(/i) (assumed constant) and the impurity level (partial) widths 
r„ = 7rt^A^(/i) for the superconductor n = 1,2,3, the basic input quantities are defined below, [with lu = ujk ^ 
{2k + l)nT (fc = 0, 1, 2, . . . ) a Matsubara frequency at temperature T = 1//3 and J2u /(^) = X^fc fi^k)]- 

3 

a„=i,2(w) = Tn[Al + + 2A„J™cos(6'i - 6I2) + w^]"^ as ^ TsiAl+uj'Y^, vi^) = ^[1 + ^cin{uj)]. 

71 = 1 

gi(ti;) = Q!i(a;)Ai + a2(w) Ji, g2(^) = a2(t^)A2 + Q!i(cj) J2, 53(^1;) = 03(0;) A3. 
3 

n— 1 7i^n' 

The self consistent equation for the field 7 (analogous to Eq. 7 in Ref?^ or Eq. 12 in RefJ^) reads, 

2U ^ [72-77(^)2 -£2- F(w,6>)]2 + 472?y(tj)2 
whose solution 7(0) is used below. 

Analysis of the results: The free energy associated with the impurity is the coefficient of — /3 in the exponent of 
the partition function Z = e^^'^, and the formalism described above yields for it the following expression, 

^^^^ = irj] ln/34{[7^ - viu;? - - F{u;, 9)]^ + ^I'viu^f}. (7) 

It will be examined as function of J and U fixing other parameters as, 

Ai,3 = 1.0, A2 = 0.8, Ti = 0.5, Ta = r3 = 0.2, £ = -2, T = 0.0005, (8) 

(energies are expressed in unit of Ai = 1). Note the ratio (5 = = 2.5 and recall that for noticeable different 
Josephson tunneling strengths Fi 7^ F2, there is no TRVGS in a 2B — IB junctioniS . 

Employing gauge invariance and setting 6*3 = 0, the minimum 6 = {61, 62, 0) of T{0) is located. If sin^i = sin ^2 = 
one evidently has a TRCGS where the Josephson current vanishes. On the other hand, when the above condition is 
not satisfied, one has a TRVGS. In the absence of Anderson impurity and close to perfect symmetry (5 « 1 it occurs 
for J > at ^1 w —62 ~ f . For U ^ the property 9i ~ —62 « § is expected to be somewhat modified, and regions 
with sin 61 , sin ^2/0 are those where TRVGS is realized. 

The phase diagram displaying these two domains is shown in Fig.[2ti- The order parameter phases of the TRVGS just 
above the border line Jc{U,5) are displayed in Fig.l^b, and indeed, they are neither nor tt so that sin^i,sin0i 0. 
The scenario emerging from Fig.[2]is as follows: For a given asymmetry parameter (5 3> 1 there is a domain in the U — J 
plane for which TRVGS is realized. It is bounded on the left and right by vertical lines U — Ui{S) and U ~ 1/2(6) and 
below by a border line Jc{U, 6) > 0. It is checked that Ui{5) < C/q = 2|e| is a slowly increasing function of 6 starting 
at Ui{l) = while U2{S) > Uq is very weakly dependent on 6. Moreover, Jc{U,6) is an increasing function of U and 
6 (actually it saturates around Uq). 

To understand the mechanism which drives the formation of the TRVGS state, we first consider the case of single- 
band superconductor Josephson junctioni^. It was found that the effective Josephson coupling between the two 
superconductors mediated through the Anderson Impurity changes sign when U goes through a critical value Uc- The 
two superconductors are in phase {S — 0) a.t U < Uc, and arc out of phase {6 = tt) at U > Uc- Associated with this 
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FIG. 2: (a) Phase diagram showing the line Jc{U,5) separating TRCGS from TRVGS in the U ~ J plane. All parameters 
(except U which varies) are defined in Eq. [S] Note the sharp cutoff points at Ui(5) ~ 3.5 (below the point B), and U2(5) ~ 6.1 
(below the point C.) (b) The phases 61,62 of the order parameters at the ground state configuration (minimum point of the 
free energy, Eq. ([7])), are displayed along the line ABCD of Fig. [2^. Along the border line Jc{U,S) we have sin ^1, sin ^1 7^ 0. 

transition is the formation of a large magnetic moment at the Anderson impurity at U > Uc- We now replace one 
of the superconductor by a two-band superconductor. In this case the effective Josephson coupling i'jos between the 
different 0„'s generated by the impurity has the form 

i^jos TisiU) cos(0i - 03) + T23{U) cos(02 - ^3) + Ti2{U) cos(0i - 62) + O(cos(0„) cos{e.m)) + 

where Ti^lU) and T23{U) become small around U ^ Uc and higher order terms in i^jos becomes important in 
determining the phase structure. The higher order Josephson terms arises from a„(a;)'s and F{u!] 0) in Eq. ([S]) which 
has much more complicated structure than the case of one-band superconductor considered in Refii^. Apparently, the 
coupling between the superconductors and the formation of tt junction for large U introduces additional frustration 
among the phases On which gives rise to TRVGS when the first order coupling terms Tia, T23 become small, i.e. around 
the critical regime U ^ Uq- 

The question now is whether this phase where TRVGS occurs can be traced experimentally A promising direction is 
to inspect the impurity magnetization m and its behavior as the phase boundary is crossed. Specifically, we fix the 
repulsive potential U and inspect the magnetization as J is varied, thereby crossing the phase boundary at Jc{U,S). 
In Ref it was found that for a IB-U-IB junction, the ground state energy and the magnetization develop kinks near 
the critical point (when studied as function of the phase difference between left and right superconductors). Here, 
on the other hand, the behavior of the ground state energy (as function of J) appears to be very smooth and rather 
slowly varying. It should be stressed that here we do not vary the phases of the order parameters, as they arc fixed 
at the ground state position {61,62). The mean field magnetization m ~ and the ground state order-parameter 
phases (61,62) are displayed as function of J for L/ = 4 in Fig.[3l The magnetization in the TRVGS phase is smaller 
than in the TRCGS phase because according to the spin analog developed in Rei^ the former phase is characterized 
by frustration. The fact that 61^2 also shows step behavior as J, and that these steps arc commensurate with those 
of the magnetization lead us to conjecture that the system has many frustrated states close in energy. 
Conclusions: The study of 2B ~U — IB carried out here is motivated by the renewed interest in multi-band super- 
conductors (stemming from the analysis of the Iron based superconductors) . The pertinent physics is fundamentally 
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FIG. 3; Impurity magnetization m (upper line) at the ground state is displayed as function of J for f/ = 4 with other 
parameters defined in Eq. (jSj. Also displayed are the phases {61,62) (divided by 2tt) of the order parameters indicating the 
ground state configuration. Both the magnetization and the phases undergo a remarkable change of behavior at the border 
line where J — Jc{U,S) (according to Fig. [2^ Jc{U = 4,5) ~ 0.1). The magnetization becomes smaller due to frustration in 
the TRVGS phase and the commensurability of steps indicates that the system passes through higher and higher frustrated 
ground-states. At the same time, it is verified that the ground state energy is a very smooth function of both J and U. 



distinct from that of a IB — U — IB systeroi^ii^ and the 2B — IB junction discussed in Reff^^ since the roles of 
strong correlations U and the coupling J between the two order parameters in the two-band superconductor interlace. 
It has been shown that for J > Jc{U,S) > and for Ui < U < U2 & TRVGS emerges which supports a non-zero 
Josephson current even without a magnetic field. The role of strong correlations as controlling the TR symmetry of 
the ground-state is evidently remarkable. As far as an experimental detection is concerned, beside the experiments 
proposed in Ref J^, our preliminary results indicate that some quantities (such as impurity magnetization) undergo a 
dramatic change as J passe through the phase boundary Jc{U, 6), hence a Josephson junction with Anderson impurity 
can serve as a potential tool for probing the relative phase of the two order parameters in a two-band superconductor 
(which is a very elusive quantity). 
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